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Abstract 



Supersymmetric grand unified models based on 50(10) gauge symmetry 
have many desireable features, including the unification of Yukawa couplings. 
Including D-term contributions to scalar masses arising from the breakdown 
of 50(10), Yukawa coupling unification only to 30% can be achieved in models 
with a positive superpotential Higgs mass. The superparticle mass spectrum 
is highly constrained, and yields relatively light top squarks and charginos. 
Surprisingly, the pattern of GUT scale soft SUSY breaking masses are close to 
those found in the context of inverted hierarchy models. Our analysis supports 
the idea that the low energy MSSM parameter space is an approximate 5O(10) 
inspired fixed point. 

PACS numbers: 14.80.Ly, 13.85. Qk, 11.30.Pb 
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Supersymmetric grand unified theories (SUSYGUTS) based on the gauge group 50(10) 
are especially attractive Not only do they unify the three forces of the Standard Model 
(SM), but they unify the matter content of each generation into a single 16 dimensional 
irreducible representation of SO(10): ip(16). The includes not only the matter su- 

perfields of the minimal supersymmetric standard model (MSSM), but also a gauge singlet 
superfield N c which includes a right handed neutrino. The gauge singlet superpotential mass 
term Mjy can be of order Mqut, and leads to sub-eV scale masses for left handed neutrinos 
(in accord with SuperK results on atmospheric neutrino oscillations), while the right handed 
neutrinos decouple via the well-known see-saw mechanism ||. In the simplest models, the 
two Higgs doublet superfields of the MSSM reside in a single ten dimensional representation 
0(10). Then 50(10) SUSYGUT models contain a superpotential interaction term 

f3fU<P+--- (1) 

where / is the Yukawa coupling which leads to masses for quarks and leptons (at this stage 
we neglect intergenerational mixing effects). Thus, the Yukawa couplings of each generation 
are assumed to be unified above the GUT scale. More sophisticated treatments of the 
Yukawa matrices can allow for predictions of all SM fermion masses and mixing angles in 
terms of just a few parameters M. Here we will focus only on third generation Yukawa 
couplings, since they will be large, and can have a substantial impact on the spectrum of 
superpartners. 

In much the same way that the three gauge couplings of the MSSM can be extrapolated 
from their weak scales values to their GUT scale values using renormalization group (RG) 
evolution, so too can the Yukawa couplings be evolved from the weak to the GUT scale to 
test models with Yukawa coupling unification. In the bottom-up approach used in ISAJET 
v7.58 we begin with the weak scale values of mj, m t and m T in the DR regularization 
scheme: m T (M z ) = 1.7463 GeV, m b (M z ) = 2.92 GeV, and m t {m t ) = 163.4 GeV. We 
calculate the corresponding Yukawa couplings, and evolve both Yukawa and gauge couplings 
to higher energies using two- loop RG equations (RGEs) ||. Once M GUT is determined (by 
the point at which the 577(2) and £7(1) gauge couplings meet), we evolve gauge couplings, 
Yukawa couplings and soft SUSY breaking (SSB) mass terms to the weak scale M wea k, where 
electroweak symmetry is broken radiatively (REWSB), and where the entire SUSY particle 
mass spectrum can be calculated. At this stage, the Yukawa couplings can be updated 
to include SM and MSSM loop corrections and the RGE process is iterated until a 

convergent spectrum of SUSY particle masses is obtained. In this way, the GUT scale values 
of the Yukawa couplings depend on the SUSY particle mass spectrum. 

In previous reports f8HT0| , unification of Yukawa couplings was investigated within the 
context of 50(10) SUSYGUT models. It is well known that in the mSUGRA model, with 
universal SSB masses at the GUT scale, a high degree of Yukawa coupling unification only 
occurs for negative values of the superpotential Higgs mass term /z, and for values of the 
ratio of Higgs field vevs tan/5 = ^ ~ 50 ||7|U- For such high values of tan/3, the scalar 
potential no longer has the appropriate form at the weak scale to accommodate REWSB. 
Generally, the SSB down Higgs mass squared gets driven to negative values before 
the up Higgs mass squared m? H . However, even if scalar masses are universal above the 
GUT scale, they will receive JJ-term mass contributions at the GUT scale arising from the 
breakdown of 5O(10) gauge symmetry. Thus, scalar masses are shifted by an amount 
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m Q = m E = m U = m 16 + M D 

m 2 D = m\ = ra 2 6 — 3Mp 

m l = m w + 5M i> 
= m ?o T 2M 2 

where M|, parametrizes our ignorance of the exact breakdown mechanism for 5*0(10), and 
can have either positive or negative values of order M 2 eak . Thus, the model is characterized 
by the following free parameters: 

m w, m io, M D , mi/ 2 , A , tan/? and sign(fi). (2) 

The value of tan j3 will be restricted by the requirement of Yukawa coupling unification to 
be close to ~ 50. In this model, for positive values of M|,, the GUT scale values of rn 2 Hu 
and m\ are split, and m 2 Hu gets a head start on running towards negative values. For a 
sufficiently large value of M|,, REWSB can be recovered, and viable supersymmetric models 
with a high degree of third generation Yukawa coupling can be generated ||. 

In previous works, model parameter space was mapped out under the restriction of GUT 
scale Yukawa unification to 5% H , and implications of the model for cosmological neutralino 
relic density, direct detection of dark matter, radiative decays b — > s'j and collider searches 
were determined [pi]]. A favorable relic density was obtained over much of model parameter 



space owing to s-channel neutralino pair annihilation via the very wide A and H Higgs poles 



at large tan (3 |Tl| . However, the decay width for b — > s-y was found to be very large and 
generally in disagreement with experimental limits unless one entered the decoupling regime, 
where model parameters began becoming unnatural. The large b — > 57 branching fraction 



at large tan/3 and \i < is well known JTSJ. However, for /i > 0, the b — > 57 branching ratio 



can be in accord with experimental limits at large tan (3. 

In addition, the recent measurement by the E821 experiment |13| of the muon anomalous 
magnetic moment a M = (g — 2)/2 was found to deviate from SM predictions [1^] by 2.6a. If 
the deviation is interpreted in terms of supersymmetric models, then it disfavors models with 



/i < 0, assuming positive SSB gaugino masses |lq| . For these reasons, it seemed prudent to 
re-examine Yukawa unification for positive values of /1, relaxing the ad-hoc 5% criteria used 
in Refs. @,[IpJ. 

Our procedure is as follows. We generate random samples of model parameters 

< m 16 < 2000 GeV, 
< mio < 2500 GeV, 
< mi/a < 1000 GeV, 
-1000 2 < M 2 D < +1000 2 GeV 2 , 

45 < tan/3 < 55, 
-5000 <A < 5000 GeV and 
/I > . 

We then calculate the non-universal scalar masses according to formulae given above, and 
enter the parameters into the computer program ISASUGRA. ISASUGRA is a part of the 
ISAJET package H which calculates an iterative solution to the 26 coupled RGEs of the 
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MSSM. We required unification of third generation Yukawa couplings at the GUT scale to 
45%. Our requirement of 45% unification is determined by defining the variables ?v, r t b and 
r tT , where for instance r hr = max (fb/ f T , fr /fb)- We then require R = max(r bT ,r tb ,r tT ) < 
1.45. We were able to find many solutions fulfilling the above criteria, with the best overall 
unification achieved for R = 1.28, or Yukawa coupling unification to 28%. The top and tau 
Yukawa couplings could frequently unify to very high precision, while fb/ ft and f b / f T could 
unify to ~ 0.7 — 0.8. The Yukawa couplings at the GUT scale, fi(McuT), should differ from 
the unified Yukawa, four, due to threshold corrections, fi(M GUT ) = /gi/t(1 + q)- It would 
be difficult to explain our 30 — 40% deviation from perfect Yukawa unification by GUT 
scale threshold corrections in simple SUSYGUT models ||16|| , where threshold corrections 
are expected of a few percent. More complicated models may be needed to explain the large 
threshold corrections needed for SO(10) SUSYGUT models with /i > |i~7j . 



The parameter space regions with Yukawa coupling unification to 45% are shown in 
Fig. [I]. Dots represent allowed solutions, while crosses denote solutions in violation of the 
experimental constraints m~± < 100 GeV, m^ < 113 GeV and m fl < 73 GeV. The values of 
tan p generated were within the narrow range of tan f3 ~ 46 — 48, typically somewhat lower 
than results assuming fi < 0. In frame a), we show models in the m\§ vs. m\/2 plane. The 
values of m 16 and m 10 typically lie beyond 1 TeV and cluster about the line m 10 = -\/2^i6- 

In frame b), we show solutions in the m 16 vs. plane. In this case, we find 
restricted to values of —0.4 to +0.4 TeV. Solutions can be obtained with Mo ~ 0, which 
brings us back to the mSUGRA model. This contrasts with the fi < case where 5% 
unification required strictly positive D-terms. 

In frame c), we show the m\Q vs. Aq parameter plane. Surprisingly, the best Yukawa 
unified solutions are found only for large negative values of the A Q parameter, and cluster 
about the line A = —2miQ. Finally, in frame d), we show the m\Q vs. m 1 / 2 plane, and find 
model solutions occuring for ~ 0.1 — 0.9 TeV, with m 16 always greater than m 1( / 2 . 

It is particularly intriguing that model solutions with the best Yukawa coupling unifica- 
tion cluster about model parameters with 

Al = 2ml = Am\,. (3) 

These particular boundary conditions were found by Bagger et al. from a rather different 



approach [IS|, by looking analytically for fixed point behavior in third generation SSB scalar 
masses which would give rise to SUSY models with a radiatively driven inverted scalar mass 
hierarchy (RIMH). In these models, one may begin with GUT scale scalar masses beyond 
the TeV scale. RG evolution drives third generation and Higgs scalar masses towards zero, 
while scalars of the first two generations remain beyond the TeV range. In this way, multi- 
TeV first and second generation scalar masses act to suppress SUSY flavor and CP violating 
processes, while still fulfilling conditions of naturalness, which mainly apply to the sub- TeV 
third generation and Higgs scalar SSB masses. The RIMH mechanism is viable for SO (10) 
based models with Yukawa coupling unification upon implementation of the above specific 



set of scalar mass boundary conditions [19|,[2(J . Our results here are obtained using a bottom- 
up approach, and indicate that for \i > 0, a high degree of Yukawa coupling unification can 
only be obtained using approximately the boundary conditions Eq. [| The values of the 
Yukawa couplings obtained at the GUT scale are f t : 0.47 — 0.50, : 0.35 — 0.37 and 



f T : 0.47 — 0.52. Their magnitudes are only sufficient to generate a small scalar IMH [ 19, 20 
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For Yukawa unified solutions with // < 0, we found a weaker correlation for the soft masses 
around mio = v2mi6 @ but we find no correlation for the trilinear parameter. 

In Fig. we show values of selected weak scale sparticle masses generated from the 
Yukawa unified model. In frame a), we show the tjia vs. rrih plane of Higgs masses. The 
light scalar h has a mass clustering about the region ~ 115 — 130 GeV, while rriA ranges 
between 100 — 500 GeV, and is generally lower than values generated in models with 
lower values of tan/3. Both the h and A (and also the heavy Higgs H) may be accessible 
to Higgs searches at the Fermilab Tevatron [plJI , and the low values of rriA may give rise to 
measureable rates for B — > ji + pT decay p2f . 

In frame b), we show the vs. mr plane. We always find m~ hx > m^, in contrast to 
Yukawa unified models with /i < 0. The value of ranges between 100 — 600 GeV for 
m~ bi < 1 TeV. 

Frame c) shows the m~± vs. m?, plane. We note that m f , > m~±, while m~± < 400 
GeV, and is likely accessible to a linear e + e~ collider operating with ~ 800 GeV p3"| , |24"|| . 

Finally, in frame cf), we show the mg R vs. rriu R plane. In this case, rrie R and m Sfl typically 
lie beyond 1 TeV. Such high mass values for first and second generation scalars can act to 
suppress many CP violating processes via a decoupling solution; they are generally not 
sufficiently heavy to suitably suppress the most dangerous flavor violating processes, such 
as K — K mixing [p5 |. A sample sparticle mass spectrum is shown in Table 1 for a Yukawa 



unified model with \i > 0. In this = 16.6 x 10~ 10 , within the 2a limit from E821 

||13|| : for our other models, typically varies around (10 ± 6) x 10~ 10 . 

The cosmological relic density of neutralinos has been calculated in Ref. 0] for Yukawa 
unified models with \i < 0. Little should change by switching to /i > 0: the pseudoscalar 
A and heavy scalar Higgs H will still have large widths of order 20 — 50 GeV due to the 
large b and r Yukawa couplings, and will be light enough that x?Xi annihilation can take 
place efficiently through s-channel annihilation. In addition, the rates for direct detection 



of relic neutralinos will remain large, as in the /i < case |I0" |. The rate for b — *> sj can 
be substantially different for /i > compared to the /i < result, and regions of parameter 
space certainly exist where this decay rate falls within experimental limits. Explicit results 



for the RIMH model with /i > have been shown in Ref. p(J . Finally, the value of a M has 



been calculated in Ref. |L5|] for Yukawa unified models with /i < and for RIMH models 
with > 0. Regions of model space with an acceptable a M certainly exist for the /i > case. 
Further results along these lines will be presented in a forthcoming publication. 

Finally, we note that Yukawa unified models with /i > have also been recently reported 



by Blazek et al. [p6] . These authors use a top-down approach and adopt independent values 



for mH u and rather than D-term splitting amongst scalar masses. We verify that in this 
case also Yukawa unified solutions can be obtained, although in our approach these generally 
occur at the 35-55% level. The solutions typically have ~ v^2^i6, with rriH u < vrtH dl 
and A ~ — 2m 16 . 

It is well known that the three standard model gauge couplings approximately unify when 
extrapolated to the scale M GUT ~ 2 x 10 16 GeV. This may be regarded as a coincidence, or as 
evidence for SUSYGUTs. Similarly, as a consequence of Yukawa unification, the clustering 
of the SSB parameters about an approximate 5*0(10) inspired fixed point, if taken seriously, 
can be considered as evidence for a supersymmetric .SO (10) grand unified model. 
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TABLES 



TABLE I. Weak scale sparticle masses and parameters (GeV) for an 50(10) case study. 



parameter value 

nm 1108.3 

mio 1497.7 

M D 58.1 

mi/2 325.2 

A -2108.0 

tan (3 49.9 

ft(M GUT ) 0.484 

MMqut) 0.372 

f T (M GUT ) 0.518 

m § 810.9 

m Sl 1267.1 

m,- 1253.0 

m fl 211.3 

m 6i 607.7 

m £ - L 1120.3 

m £ ~^ 1115.4 

m* 1117.5 

m fl 308.1 

mo T 829.1 

m~± 145.5 

m~o 164.6 

m~o 111.6 

Xi 

m h 121.6 

200.1 
225.6 

H 165.4 



16.6 x 10" 10 
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FIG. 1. Plots of regions of parameter space where valid solutions to minimal SUSY SO(W) 
are obtained, consistent with Yukawa coupling unification to 45% with fj, > and REWSB. Crosses 
are excluded by collider searches. 
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FIG. 2. The range of selected sparticle masses that are generated in minimal SUSY SO(W) 
models with Yukawa coupling unification to 45%, \x > and REWSB. Crosses are excluded by 
collider searches. 
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